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Abstract 

By first principle calculations, a dual spin filter effect under finite bias voltages is demonstrated 
in an antiferromagnetic junction of symmetric zigzag graphene nanoribbon (ZGNR). Unlike con- 
ventional spin filter devices using half metallic materials, the up- and down-spin electrons are 
unidirectionally filtered in the counter direction of the bias voltage, making the junction a dual 
spin filter. On the contrary, asymmetric ZGNRs do not exhibit such a spin filter effect. By ana- 
lyzing Wannier functions and a tight-binding model, we clarify that an interplay between the spin 
polarized band structure of tt and vr* states near the Fermi level and decoupling of the interband 
hopping of the two states, arising from the symmetry of the wave functions, plays a crucial role in 
the effect. 

PACS numbers: 85.75.-d, 72.25.-b, 73.63.-b 
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I. INTRODUCTION 



The graphenes have been recently attracting much attention as a candidate material for 
spintronics devices because of its peculiar electronic structures,-'^'^'^'^'^'^ while so far most of 
promising materials in developing the devices have been found in ferromagnetic (FM) met- 
als and FM semiconductors such as GaMnAs.— In fact, recent experiments have achieved 
spin injection into graphene layers at room temperature for the first time among molecular 
materials,->2 and observed the magnetoresistance (MR) effect.- In addition, several intrigu- 
ing transport properties have been theoretically predicted especially for zigzag graphene 
nanoribbons [ZGl^Ks) M^'^'^^^^^^^^^^ For instance, it is shown that ZGNRs can be utilized 
even for generation of a spin polarized current in which an external electric field applied 
across ZGNR along the lateral direction may induce a half-metallic band structure being 
responsible for the spin polarized current.- Also an extraordinary large MR and generation 
of spin-polarized current are also theoretically found in a spin valve device consisting of 
ZGNR.- The unique properties in the electronic transport properties of ZGNRs can be at- 
tributed to the characteristic band structures of ZGNR and the symmetry of wave functions 
near the Fermi level.-i^i^ Along this line, the family of graphenes related to ZGNR might 
be anticipated to exhibit unexpected electronic properties which can be useful in developing 
spintronics devices. However, transport properties of ZGNRs under finite source-drain bias 
voltages Vbias have not been fully explored even from the theoretical points of view. For de- 
velopment of future device applications it is highly important to fully understand intrinsic 
behaviors in the transport properties of ZGNRs under a wide range of finite bias voltages. 

In this paper we present a novel and intrinsic electronic transport property of ZGNR 
under finite bias voltage Hias, based on first principle calculations, that a symmetric ZGNR 
with an antiferromagnetic (AFM) junction exhibits a dual spin filter effect under finite 
Vbias; namely, the up- and down-spin electrons are unidirectionally filtered in the counter 
direction of Vuas- Based on an analysis using Wannier functions (WFs) and a tight-binding 
(TB) model, we further clarify that the spin filter effect arises from an interplay between 
the band structure and absence of the interband hopping of the vr and vr* states near the 
Fermi level. 
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II. COMPUTATIONAL DETAILS 



At each bias voltage Vbias the electronic structure of ZGNR shown in Fig. 1(a) is self- 
consistently determined under a temperature of 300 K by means of a non-equilibrium Green 
function (NEGF) methocy^>i^ coupled with a local spin density approximation (LSDA)^^ in 
the density functional theory (DFT).— The equilibrium density matrix (DM) is evaluated 
by a contour integration method/^ while the non-equilibrium DM is numerically computed 
on a line with an imaginary part of 0.01 eV which is parallel to real axis in complex plane.— 
The conductance and current are calculated by the Landauer formula.— Norm-conserving 
pseudopotentials are used in a separable form with multiple projectors to replace the deep 
core potential into a shallow potential.— Pseudo-atomic orbitals (PAOs) centered on atomic 
sites are used as basis functions. The PAO basis functions, generated by a confinement 
scheme,— are specified by H5.5-s2 and C4.5-s2p2, where the abbreviation of basis functions, 
such as C4.5-s2p2, represents that C stands for the atomic symbol, 4.5 the cutoff radius 
(Bohr) in the generation by the confinement scheme,— and s2p2 means the employment of 
two primitive orbital for each of s- and p-orbitals. The real space grid techniques are used 
with the energy cutoff of 120 Ry as a required cutoff energy in numerical integrations and 
the solution of Poisson equation using FFT.— In addition, the projector expansion method is 
employed in the calculation of three-center integrals for the deep neutral atom potentials.— 
The geometrical structures used are optimized with a criterion of 10~^ hartree/bohr for 
forces on atoms under the periodic boundary condition. All the calculations were performed 
by an ab initio DFT code, OpenMX.— 



III. RESULTS 



77,-ZGNRs are characterized by the number of carbon atoms, even (symmetric) or odd 
(antisymmetric), in the sublattice being across ZGNR along the lateral direction, while 
Fig. 1(a) shows the case of eight carbon atoms abbreviated as 8-ZGNR. For ZGNRs we focus 
a spin configuration with an AFM junction as shown in Fig. 1(a). Considerable magnetic 
moments are found at both the zigzag edges which can be attributed to the existence of the 
fiat band near X-point.— Although the AFM coupling between the zigzag edges is favored 
by about 10 meV per edge carbon atom compared to the FM coupling, we consider the 
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spin configuration consisting of tlie FM coupling between tlie zigzag edges and tlie AFM 
junction at tlie central region. The spin configuration is crucial for the spin filter effect 
we discuss, and might be realized by a magnetic field applied in a spin valve device or 
chemical modifications for ZGNRs. It should be noted that the effect does not appear in 
the other spin configurations. The current-voltage (I-V) characteristics for 8-ZGNR with 
the spin configuration depends on intriguingly not only spin, but also the direction of Vbias 
as shown in Figs. 1(b) and (c). Interestingly, the up-spin electron substantially flows only in 
the negative regime of Vbias, while the down-spin electron flows only in the positive regime. 
The ratio of the spin dependent currents, /up//down, turns out to be 44.3 at —0.4 V, and 
the value is equivalent to the rectification ratio for each spin-dependent current because of 
the I-V characteristics. The I-V characteristics of other even cases, n =6 and 10, are also 
found to be nearly equivalent. Thus, a symmetric ZGNR plays dual roles as a unidirectional 
spin filter for each spin state under finite Vbias- The effect can also be regarded as a dual 
spin diode effect due to the unidirectional nature of the spin dependent current. In contrast, 
the current for 7-ZGNR with the spin configuration is almost independent of spin, and 
proportional to Vbias within the regime, leading to /up/ /down of a nearly one. Thus, it should 
be emphasized that the parity in the geometrical structure of n-ZGNR is the key factor even 
for the behavior of the spin polarized current as in the spin unpolarized current .- 

The I-V characteristics for 8-ZGNR is understood by examining the dependency of the 
conductance on Vbias- It is found in Fig. 2(a) that there is a conductance gap of about 0.4 
eV for both the up- and down-spin states around the chemical potential at Vbias = V. 
The gap for the up-spin state decreases as Vbias increases toward the negative direction, and 
approaches zero at about Vbias = —0.4 V as shown in Fig. 2(b), while the zero gap is kept up 
to Vbias = —0.55 V. From then onwards the gap increases as illustrated in Fig. 2(c). On the 
other hand, the gap for the down-spin state monotonically increases as Vbias increases toward 
the negative direction. The conductance gaps for the up- and down-spin states are shown as 
a function of Vbias in the inset of Fig. 2(a). We see that the behavior of the gaps is opposite 
for the up- and down-spin states in the regime of the positive direction of Vbias in contrast 
with that in the negative direction. The shaded regions in Fig. 2(b) and (c) correspond to 
the regime in between two chemical potentials of the left and right leads, where the electron 
transmission can contribute to the current flow. Thus, it is confirmed that in the regime 
of the negative Vbias, l-^upl increases up to around Vbias = —0.55 V, where the gap opening 
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occurs again, and from then onwards saturated. Also it turns out that |/down| should be 
nearly zero in the regime of the negative Vbias, since the shaded region is always inlying in 
the conductance gap for the down-spin state. In the regime of the positive Vbias, the I-V 
characteristics can be explained in the same way by considering the reversed roles of up- and 
down-spin states. Note that the reversal of the roles arises only in the spin configuration 
with the AFM junction. 

The opening and closing of the conductance gap can be attributed to the band structure 
near the Fermi level of 8-ZGNR with the FM coupling between the zigzag edges. The 
exchange splitting A^; is found to be 0.553 eV at X-point from Fig. 3(a). The right blue and 
pink shades denote the conductance gaps for the up- and down-spin states in Fig. 2. It is 
clearly seen that the conductance gap corresponds to an energy regime where the 7r(7r*) state 
in the left panel overlaps with only the 7r*(7r) state with the same spin in the corresponding 
right panel. The vr* up-spin state in Fig. 3(a) overlaps with only the vr up-spin state in 
Fig. 3(d) in the conductance gap. Once the overlap regime fades away at the situation given 
by Figs. 3(b) and (e), and then onward it turns out that the vr up-spin state in Fig. 3(c) 
overlaps with only the vr* up-spin state in Fig. 3(f) in the regime of the conductance gap. 
On the other hand, for the down-spin state the energy regime, where the tt state in the left 
panel overlaps with only the vr* state in the right panel, linearly increases as the energy shift 
increases toward the negative direction as shown in Figs. 3(d), (e), and (f). The same idea 
except for reversing of the roles of the vr and vr* states can apply to the case of the energy 
shift toward the positive direction. Since the band structures in the left and right panels can 
be regarded as those of the left and right leads in Fig. 1(a), the correspondence implies that 
the electron transmission from the 7r(7r*) to tt*{tt) states is forbidden. In fact, it is shown 
that the Bloch function of the vr (vr*) state is antisymmetric (symmetric) with respect to the 
a mirror plane which is the mid plane between two edges.-i^ Thus, the electron transmission 
should be forbidden in the energy regime where the vr (vr*) overlaps with only the vr* (vr) 
states. The band structure of 7-ZGNR with the FM coupling between the zigzag edges 
is very similar to that of 8-ZGNR so that the above analysis can also apply to the case. 
However, as discussed later the transmission is allowed even for the energy regime, where 
only the overlap between the vr (vr*) and the vr* (vr) states survives, due to the absence of 
the mirror plane in 7-ZGNR, leading to the linear I-V characteristics as shown in Figs. 1(b) 
and (c). 
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To further verify the physical origin, explained above, of the peculiar I-V curves which 
can be determined by the characteristic band structure of ZGNRs and the symmetry of 
wave functions, we construct WFs for the vr and vr* states for 7- and 8-ZGNRs with the non- 
magnetic state,^ and evaluate tight-binding (TB) parameters using WFs. It is confirmed 
that WFs of 7-ZGNR is neither symmetric nor antisymmetric (not shown), and that the 
absolute value of the nearest neighbor hopping integral between the vr and vr* states is 
0.485 eV which is comparable to the nearest neighbor hopping integral for the vr (vr*) state 
of -0.777 (0.784) eV. On the other hand, WFs for the vr and vr* states of 8-ZGNR are 
antisymmetric and symmetric with respect to the a mirror plane, and localized in three unit 
cells along the ribbon direction as shown in Figs. 4(a) and (b). We confirm that the hopping 
integrals between WFs for the vr and tt* states are nearly zero (not shown). Also it turns 
out that the next nearest neighbor hopping integral and more, /i„(n = 2 — 4), are one order 
smaller than hi as shown in Table I. Thus, we can construct a simple TB model determined 
by a fitting, containing the essence of the spin filter effect, which consists of only the on-site 
energy, the nearest neighbor hopping integral, and the exchange splitting A^. The fitted 
parameters in Table I are determined so that the band width and the degeneracy of the 
TT and vr* states at X-point can be reproduced. It should be noted that the two band TB 
model can be decomposed into two decoupled one band models due to the decoupling of the 
TT and 71* states. The fact allows us to evaluate a spin polarized current for the 7r(7r*) state 
using the Landauer formulai^i for the one band TB model and the fitted parameters as: 



dE{fL-fR)T{E), 



(1) 



where fi and are the Fermi functions with the chemical potentials for the left and right 
leads, respectively. The transmission T is given by an analytic formula: 

ASL{E)Sa{E) 



T{E) 



[SL{E) + Sn{E)y 



with 



Sl{E) 




(2) 



(3) 



(4) 



where Sl and Sr are not zero only if the number in the square root is larger than zero, and 
therefore the transmission is finite within the energy range where both Sl and Sr survive. It 



is emphasized that the behaviors of the conductance gap in Fig. (2) can be easily traced by 
the model. Using Eq. (1) we evaluate the gate voltage dependency of the I-V characteristics 
for the up-spin state as shown in Fig. 4(c), where A3, of 0.553 eV, being that of the spin 
polarized 8-ZGNR, is used, and the gate voltage legate is taken into account by adding legate 
to the on-site energy. The proposed TB model can be validated by the fact that the I-V 
characteristics at V^ate = 0.0 V is quite similar to that in Fig. 1(b) calculated by the NEGF 
method. In addition, it is found that the I-V characteristics does not largely change as long 
as Vgaxe lies in between i^A^., while the absolute threshold voltage at which the current 
starts to flow can be tuned. However, once | Viatel exceeds ^A^,, the originally suppressed 
current starts to leak in the positive Vbias regime. The leak current is due to the fact that 
one of the chemical potentials is located at the outside of the conductance gap. In case 
that ZGNR contacts with a metallic substrate, the Fermi level of ZGNR might be shifted 
by charge transfer between them. For such a case the gate voltage dependency apparently 
suggests that the Fermi levels must be adjusted within the energy regime of the exchange 
splitting by applying a proper gate voltage in order to keep the high /up //down- 



IV. CONCLUSIONS 



In summary, we demonstrate based on the NEGF method coupled with DFT that the 
symmetric ZGNRs with an AFM junction possess intrinsically a peculiar I-V characteristics 
which can be regarded as a dual spin filter effect. It is shown by analyzing the band structure, 
WFs of the TT and tt* states, and a TB model that the physical origin of the spin filter effect 
can be attributed to the spin polarized band structure of the symmetric ZGNR and the 
symmetry of wave functions of the tt and vr* states near the Fermi level. The dual spin filter 
effect might initiate a novel avenue in developing spintronics using graphene based devices. 
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Table: 



TABLE I: Tight-binding parameters (eV) evaluated by WFs denoted by WF, and a fitting, 
denoted by fitted, for the tt and tt* states of the non-spin polarized 8-ZGNR, where £ is the on-site 
energy, and hi, h2 - ■ ■ the nearest and the second nearest neighbor hopping integrals, and so on. 
The Fermi level is set to zero. 

£ hi /i2 ^3 ^4 

TT (WF) -1.3609 -0.7660 0.0076 0.0529 -0.0352 
TT* (WF) 1.4486 0.7708 -0.0400 -0.0513 0.0269 

TT (fitted) -1.4165 -0.7083 

TT* (fitted) 1.4135 0.7067 
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Figure captions: 



Fig. 1 

(a) 8-ZGNR, terminated by hydrogen atoms, with an AFM junction together with the 
spatial distribution of the spin density at Vbias = V, where the red and blue colors stands 
for positive and negative signs, and L, R, and C means the left and right leads, and the 
central scattering region. /-H>ias curves for (b) the up-spin and (c) the down-spin states of 

7- and 8-ZGNRs with the AFM junction. 

Fig. 2 

Conductance of 8-ZGNR with the AFM junction with Vbias of (a) 0.0, (b) -0.4, and (c) 
—1.0 V. The inset in (a) shows the conductance gap near the chemical potentials. 

Fig. 3 

(a)-(c) Band structures of 8-ZGNR with a FM couphng corresponding to the left lead 
region in Fig. 1(a). Band structures shifted with (d) 0.0, (e) —0.4, and (f) —1.0 eV of 

8- ZGNR with a FM coupling corresponding to the right lead region in Fig. 1(a). The blue 
and red shades show the conductance gap for the up- and down-spin electrons, respectively. 
The horizontal dot line indicates the chemical potential, being position dependent, of 
8-ZGNR shown in Fig. 1(a). 

Fig. 4 

Wannier functions for (a) the tt and (b) tt* states of the non-spin polarized 8-ZGNR. 
(c) /-Hias curves of the up-spin state calculated by the simple TB model with the AFM 
junction for a series of the gate voltage V^ate- 
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